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Electromagnetic  wave  structures  within  subauroral  polarization 
streams 

E.  V.  Mishin,*  W.  J.  Burke,^  C.  Y.  Huang,*  and  F.  J.  Rich^ 

Received  25  November  2002;  revised  25  March  2003;  accepted  1  May  2003;  published  1  August  2003. 

[i]  We  report  on  oscillations  in  electric  (bEy)  and  magnetic  (6Bz)  fields  and  plasma 
density  (bNi)  observed  by  Defense  Meteorological  Satellite  Program  (DMSP)  satellites 
within  fast  subauroral  convection  streams  in  the  evening  sector  during  the  magnetic  storm 
of  6  November  2001.  There  are  two  types  of  wave  phenomena.  The  first  and  more 
common  is  characterized  by  electromagnetic  and  plasma  density  variations  that  have  the 
same  firequency  range  of  <^0.15  Hz  in  the  spacecraft  frame  of  reference.  The  second  is 
characterized  by  large-amplitude  plasma  and  field  oscillations  over  a  broader  range  of 
frequencies  ~0. 1  to  0.3  Hz.  In  this  case  the  perturbation  densities  and  fields  appear  to  have 
different  frequency  responses.  In  this  and  other  magnetic  storms,  strong  waves  are 
associated  with  the  precipitation  of  ~30  keV  ions.  Ratios  of  hEylhBz  indicate  encounters 
with  mixtures  of  electromagnetic  (in  part  Alfvenic)  and  electrostatic  modes.  Poynting 
vectors  associated  with  the  oscillations  can  be  directed  either  into  or  out  of  the  ionosphere. 

The  density  perturbations  appear  to  be  extended  east-west  corrugations  in  the  plasma  flow 
streams  with  norfli-south  wavelengths  of  ~50  km.  The  8£’yand  hNt  variations  were 
anticorrelated,  as  required  for  current  conservation.  Our  analysis  shows  that  Alfvenic 
perturbations  are  consistent  with  expected  effects  of  irregular  potential  distribution  around 
ionospheric  density  irregularities  mapped  to  the  magnetosphere.  Inertial  currents  act  to 
generate  mesoscale  field-aligned  currents  carried  by  Alfven  waves,  as  was  previously 
discussed  with  regards  to  auroral  arcs  formation.  We  suggest  that  6A^  irregularities 
observed  by  DMSP  satellites  in  the  evening  sector  began  as  striated  plasma  patches  in  the 
polar  cap  that  convected  to  subauroral  latitudes.  Index  Terms:  2431  Ionosphere:  ionosphere/ 

magnetosphere  interactions  (2736);  2439  Ionosphere:  Ionospheric  irregularities;  2471  Ionosphere:  Plasma 
waves  and  instabilities;  2716  Magnetospheric  Physics:  Energetic  particles,  precipitating;  2778  Magnetospheric 
Physics:  Ring  current 

Citation;  Mishin,  E.  V.,  W.  L  Burke,  C.  Y.  Huang,  and  F.  J.  Rich,  Electromagnetic  wave  structures  within  subauroral  polarization 
streams,  /  Geophys,  Res.,  ]08(ASX  1309,  doi:  10.1 029/2002 JA009793,  2003. 


1.  Introduction 

[2]  Subauroral  magnetic  field  lines  map  to  the  inner 
magnetosphere  earthward  of  the  electron  plasma  sheet. 
During  intensifications  of  geomagnetic  activity,  electric 
fields  permeate  the  inner  magnetosphere  with  widespread 
consequences  for  the  global  distributions  of  plasmas,  par¬ 
ticles,  and  fields.  The  domain  of  auroral  particle  precipitation 
in  the  ionosphere  expands  equatorward.  In  large  numbers, 
ring  current  (RC)  ions  enter  previously  forbidden  parts  of  the 
plasmasphere.  The  source  of  the  ions  can  be  either  substorm 
injections  or  large  convection  fields  associated  with  mag¬ 
netic  storms.  Ultimately,  plasma  pressure  grows  in  the  inner 
magnetosphere.  In  regions  where  gradients  in  plasma  pres¬ 
sure  become  misaligned  with  gradients  in  the  magnetic  flux 
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tub  volume,  field-aligned  currents  (FACs)  flow  to  couple  the 
magnetosphere  with  the  conjugate  ionosphere  [Vasyliunas, 
1970].  The  recently  developed  kinetic  transport  models  of 
Fok  et  al  [2001]  and  Liemohn  et  al.  [2001]  indicate  that  the 
azimuthal  pressure  gradient  and  radial  flux  tube  volume 
gradient  are  largely  responsible  for  the  region  2  FACs.  FACs 
affect  the  distributions  of  electric  potential  required  to  satisfy 
cuirent  continuity  in  the  ionosphere.  Via  field-line  mapping, 
potential  in  the  inner  magnetosphere  is  redistributed.  Altered 
potential  distributions  are  usually  described  in  terms  of 
partial  shielding  of  the  convection  electric  field  by  the  ring 
current  [e.g.,  Southwood  and  Wolf,  197S;  Harel  et  al.,  1981; 
Spiro  et  al.,  1981;  Senior  and  Blanc,  1984;  Erickson  et  al, 
1991;  Peymirat  and  Fontaine,  1994]. 

[3]  Owing  to  the  low  conductivity  of  the  nightside 
ionosphere  at  subauroral  latitudes,  channels  of  intense, 
poleward  directed  electric  fields  develop  to  produce 
enhanced  (>500  m/s)  streams  of  westward  convection.  This 
phenomenon  is  described  by  the  names  polarization  jets  (PJ) 
[Galperin  et  al,  1974]  or  subauroral  ion  drifts  (SAID) 
[Spiro  et  al,  1979;  Anderson  et  al,  1993,  2001].  Yeh  et 
al  [1991]  and  Burke  et  al  [1998,  2000]  pointed  out  that 

.i 

1-1 


SIA  1-2 


MISHIN  ET  AL.:  SUBAURORAL  WAVE  STRUCTURES 


latitudinal  extents  of  westward  drifting  plasma  at  subauroral 
latitudes  during  periods  of  high  activity  may  be  much 
greater  than  the  to  2°  typical  of  SAID/PJ  structures. 
Recently,  Foster  and  Burke  [2002]  suggested  the  name 
subauroral  Polarization  Streams  (SAPS)  to  include  plasma 
flow  events  that  have  both  broad  and  narrow  extents  in 
latitude. 

[4]  Previous  experimental  and  modeling  efforts  treated 
SAPS  electric  fields  as  varying  smoothly  across  the  struc¬ 
ture.  However,  observations  fi*om  Defense  Meteorological 
Satellite  Program  (DMSP)  satellites  and  Millstone  Hill  radar 
[Erickson  et  al,  2002]  often  indicate  the  presence  of 
irregular  substructures  with  scale  size  of  order  tens  km 
within  SAPS.  While  Maynard  et  al.  [1980]  mentioned  that 
storm  time  electric  fields  adjacent  to  the  plasmapause  were 
strongly  variable,  consequences  of  their  observations  were 
not  pursued.  On  the  other  hand,  it  is  known  that  short-scale 
density  irregularities  in  the  F  region  can  significantly 
modify  radio  wave  propagation.  For  example,  the  HAARP 
HF  heating  facility  (see,  e.g.,  [Pedersen  and  Carlson,  2001] 
or  http://www.haaip,alaska.edu)  is  located  at  normally  sub¬ 
auroral  latitudes.  The  presence  of  plasma  irregularities  in  the 
overhead  ionosphere  can  affect  the  structure  of  the  heated 
region. 

[5]  Plasma  irregularities  of  sufficient  scale  size  can  affect 
the  dynamics  of  warm  plasmas  in  the  conjugate  magneto¬ 
sphere  and  ionosphere  [Ogawa  and  Sato,  1971].  This 
subject  has  been  addressed  in  numerous  studies  of  magne¬ 
tosphere-ionosphere  coupling  at  auroral  latitudes  [Atkinson, 
1971;  Maltsev  et  al,  1977;  Miura  and  Sato,  1980].  To  the 
degree  that  plasma  density  irregularities  represent  conduc¬ 
tance  irregularities,  polarization  electric  fields  and  field- 
aligned  currents  must  develop  to  maintain  current  continuity 
across  the  irregularities.  Modeling  by  Miura  and  Sato 
[1980]  shows  that  consequently  altered  potential  patterns 
create  feedback  conditions  inducing  structured  auroral  elec¬ 
tron  precipitation  that  enhances  the  regions  of  initially  high 
ionospheric  conductance.  Data  presented  here  mostly  come 
fi'om  subauroral  latitudes  where  no  measurable  electron 
precipitation  was  present.  We  suggest  that  an  analogous 
electromagnetic  feedback  develops,  in  this  case  between 
ionospheric  irregularities  and  conjugate  ring  current  ions. 

[6]  The  purpose  of  this  report  is  to  extend  understanding 
of  irregularity  structures  found  within  SAPS  structures.  This 
is  accomplished  through  analyses  of  perturbations  in  local 
plasma  densities  and  drifts  as  well  as  magnetic  fields 
observed  by  sensors  on  the  DMSP  FI 5  satellite  during  the 
magnetic  storm  of  6  November  2001 .  The  following  section 
describes  the  relevant  sensor  suite  on  DMSP  FI 5.  The 
observation  section  has  three  parts.  A  brief  description  of 
the  magnetic  storm  provides  a  context  for  understanding 
DMSP  measurements.  A  detailed  exposition  of  measure¬ 
ments  acquired  during  one  of  the  S  APS  crossing  illustrates 
the  quasi-dc  circuit  and  embedded  small-scale  plasma  and 
fields  irregularities.  Our  analysis  demonstrates  that  the 
irregularities  often  have  narrow  frequency  spectra.  Both 
quasi-electrostatic  and  electromagnetic  wave  structures  can 
coexist  in  close  physical  proximity.  At  times  the  Poynting 
flux  carried  by  the  electromagnetic  waves  was  directed 
away  fi-om  the  ionosphere.  Characteristics  of  the  full  set 
of  DMSP  measurements  are  summarized  in  tabular  form. 
The  final  section  discusses  possible  wave-generation  mech¬ 


anisms  and  the  roles  of  ring  current,  plasmasphere,  and 
ionosphere. 

2.  Instrumentation 

[7]  DMSP  satellites  are  three-axis  stabilized  spacecrafts 
that  fly  in  circular,  sun-synchronous  polar  (inclination 
98.7°)  orbits  at  an  altitude  of  ~840  km.  The  geographic 
local  times  of  the  orbits  are  either  near  the  1800-0600 
(F13)  or  2100-0900  (F12,  14,  15)  meridians.  Owing  to  the 
offset  between  the  geographic  and  geomagnetic  poles 
DMSP  satellites  sample  wide  range  of  magnetic  local  times 
(MLT)  over  the  course  of  a  day.  The  ascending  nodes  of 
DMSP  orbits  are  on  the  duskside  of  the  Earth.  Thus  the 
satellites  move  toward  the  northwest  in  the  evening  LT 
sector.  Each  satellite  carries  a  suite  of  sensors  to  measure  (1) 
fluxes  of  auroral  particles  (SSJ4),  (2)  the  densities,  temper¬ 
atures,  and  drift  motions  of  ionospheric  ions  and  electrons 
(SSIES),  and  (3)  perturbations  of  the  Earth  magnetic  field 
(SSM). 

[s]  Identical  SSJ4  sensors  are  mounted  on  the  top  sides  of 
DMSP  satellites  to  measure  fluxes  of  precipitating  electrons 
and  ions  in  the  energy  range  between  30  eV  and  30  keV 
[Hardy  et  al,  1984].  The  measurements  are  made  by  four 
detectors,  one  high  energy  detector  and  one  low  energy 
detector  for  the  each  of  the  particle  types.  The  ion  detectors 
have  no  mass  discrimination  capabilities.  Each  detector  has 
10  logarithmically  spaced  energy  steps.  The  high  energy 
detectors  step  from  30  keV  to  1  keV  and  the  low  energy 
detectors  step  from  1  keV  to  30  eV.  The  last  step  of  the  high 
energy  detector  is  set  at  the  same  energy  as  the  first  step  of 
the  low  energy  detector.  Only  particles  within  an  energy 
band  of  approximately  10%  of  the  channel  step  energy 
freely  pass  from  aperture  to  the  channeltron  detector.  The 
particle  fluxes  are  measured  within  a  solid  angle  of  2°  x  5° 
for  the  high  energy  channels  and  4°  x  5°  for  the  low  energy 
channels  centered  on  local  vertical.  Each  detector  has  a 
dwell  time  of  0.098  s  and  a  0.002  s  period  between  steps  to 
stabilize  the  voltage.  Each  detector  make  a  complete  10  step 
sequence  in  1  s.  One  20  point  ion  and  one  20  point  electron 
is  returned  once  per  second.  Actually,  the  spectra  have  only 
19  points  because  of  the  redundancy  of  the  1  keV  energy 
step.  Unfortunately,  fluxes  of  ions  with  energies  below 
1  keV  are  unavailable  from  the  FI 5  satellite  as  the  low- 
energy  ion  detector  failed. 

[9]  SSIES  sensors  are  mounted  on  the  ram  facing  surfa¬ 
ces  of  the  satellites.  They  consist  of  (1)  an  ion  drift  meter  to 
measure  the  horizontal  {Vh)  and  vertical  (F^)  cross-track 
components  of  plasma  drift  within  the  range  of  ±3000  m/s 
and  a  one-bit  resolution  of  12  m/s  for  ambient  ion  densities 
greater  than  5-10^  cm”^;  (2)  retarding  potential  analyzer 
to  measure  ion  temperatures  (7}),  composition,  and  the 
in-track  component  of  plasma  drift  (Fjj);  (3)  an  ion  trap 
to  measure  the  total  ion  density  {N^\  and  (4)  a  spherical 
Langmuir  probe  mounted  on  an  80-cm  boom  to  measure 
the  density  {N^  and  temperature  {T^  of  ambient  electrons 
[Rich  and  Hairston,  1984].  The  drift  measurements  are 
made  at  a  time  resolution  of  1/6  s  and  the  density  measure¬ 
ments  are  at  24  Hz. 

[10]  SSM  sensors  are  triaxial  fluxgate  magnetometers  that 
are  mounted  on  the  bodies  of  the  FI 2-F 14  spacecrafts.  For 
DMSP  flight  F15  the  magnetometers  are  mounted  on  5-m 
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UT,  hour 

Figure  1.  (top)  IMF  By  (heavy  line),  Bx  (light  line),  and 
Bz-  (bottom)  SYMH  and  AE  indices;  bold  lines  designate 
the  intervals  where  SAPS  were  observed. 


booms.  These  sensors  are  much  less  susceptible  to  space- 
craft-generated  electromagnetic  contamination  than  their 
body-mounted  predecessors.  The  high-precision  magnetic 
field  measurements  required  for  the  present  study  all  come 
from  the  F 15  satellite.  Magnetic  field  components  are 
sampled  at  a  rate  of  12  (7  and  Z)  and  10  (X)  s”^  and  are 
presented  in  a  satellite-centered  coordinate  system.  The 
X  axis  points  in  the  downward  direction.  The  7  axis  points 
along  the  spacecraft  velocity.  The  Z  axis  completes  a  right- 
hand  coordinate  system.  Positive  Z  components  generally 
point  in  the  antisunward  direction.  Data  are  presented  as 
differences  (AB)  between  measured  values  and  those 
assigned  by  the  ICjRF-90  magnetic  field  model  (B©).  Nearly 
simultaneous  observations  from  the  SSM,  SSIES,  and  SSJ4 
sensors  on  the  other  spacecrafts  are  used  to  interpret  F15 
measurements. 

3.  Observations 

[ii]  Salient  features  of  the  6  November  2001  magnetic 
storm  are  summarized  in  Figure  1.  The  top  two  plots 
respectively  show  the  GSM  Y,  X^  and  Z  components  of 
the  interplanetary  magnetic  field  (IMF)  as  measured  at  the 
ACE  satellite  near  the  first  Lagrangian  point  (Ei).  For 


comparison  with  Earth-based  measurements,  ACE  data 
were  shifted  by  one  1  to  allow  for  propagation  from  Ly 
At  the  time  of  this  writing,  solar  wind  data  for  most  of  the 
day  are  unavailable.  We  note,  however,  that  between  the 
beginning  and  end  of  the  day  the  solar  wind  speed  rose  from 
^430  to  ^^730  km/s  and  the  density  decreased  from  -^10  to 
~2  cm“^.  The  bottom  two  plots  respectively  show  the 
SYMH  [KYOTO]  and  AE  indices.  Information  provided 
by  the  SYMH  index  is  similar  to  Dst  but  is  presented  at  a 
cadence  of  1-min  versus  l-hoiu*  averages. 

[12]  Data  in  Figure  1  show  that  between  0150  and  0153 
UT  SYMH  and  AE  increased  from  —62  to  25  and  '^lOOO 
to  ~3300  nT,  respectively, .  then  rapidly  decreased  to 
~  —300  and  ^500  nT.  Shortly  after  0300  UT  SYMH 
remained  near  —300  nT  before  beginning  to  recover  1  hour 
and  20  min  later,  likely,  following  a  northward  turning  of 
the  IMF  Bz-  At  the  same  time  AE  sharply  increased  to 
~2000  nT  until  --0545  UT  The  SYMH  recovery  pro¬ 
ceeded  smoothly  to  170  nT  at  1400  UT  when  SYTMH 
again  decreased  to  '^—215  nT  for  about  2  hours.  The 
main-phase  decrease  in  SYMH  corresponds  to  a  southward 
turning  of  the  IMF  Bz  that  reached  consistent  values  near 
—70  nT.  The  second  southward  turning  of  the  IMF  Bz  to 
about  -30  nT  between  1400  and  1600  UT  corresponds  to 
the  second  SYMH  intensification.  During  this  period  AE 
experienced  ±200  nT  variations  until  about  1540  UT  when 
it  increased  up  to  ~1000  nT,  then  decreased  and  increased 
again  at  about  1730  UT,  Finally,  attention  is  directed  to 
two  bold  horizontal  lines  extending  from  about  0400  to 
0800  and  1400  to  2300  UT.  They  delineate  intervals  dur¬ 
ing  the  recovery  phase  when  FI 5  crossed  SAPS  structures 
(16  events). 

[13]  Figure  2  provides  an  example  typical  of  SAPS 
structure  crossed  by  DMSP.  It  was  sampled  during  a 
Northern  Hemisphere  pass  in  the  evening  sector  in  the 
storm’s  recovery  phase.  The  top  two  panels  show  directional 
differential  number  fluxes  of  downcoming  electrons  and 
ions,  respectively.  The  trace  in  the  third  panel  gives  1-s 
averages  of  the  horizontal  component  of  the  cross-track  drift 
velocity  {Vh)  measured  by  the  drift  meter  in  meters  per 
second.  The  light  and  heavy  traces  in  the  bottom  panel  give 
1-s  averages  of  Aff^and  AB^,  respectively.  Data  are  plotted 
as  functions  of  UT  and  MLat.  Magnetic  local  time 
(2118  MET)  practically  did  not  change.  By  definition,  the 
SAPS’  poleward  boundary  coincides  with  the  equatorward 
boundary  of  auroral  electron  precipitation  (the  dashed  line) 
and  was  crossed  at  ~1802:50  UT.  The  maximum  plasma 
velocity  of  1.75  km/s  was  sampled  at  ^1802:30  UT 
(^52°  MLat).  Irregular  structures  appear  in  and  mag¬ 
netic  field  traces. 

[14]  Figure  3  provides  the  in-track  (T)  component  of  the 
electric  field  E'y[mV/m]  =  10”VH[ni/s]  •  Box[nT],  the 
cross-track  component  of  the  magnetic  field  /XBz,  and  ion 
density  Nf.  At  northern  latitudes,  positive  slopes  in  the 
trace  correspond  to  field-aligned  currents  into  the  iono¬ 
sphere.  Apparently,  the  maximum  electric  field  roughly 
coincides  with  the  deepest  depletion  in  the  Nt  trace.  All 
three  quantities  contain  small-scale  irregularities  superposed 
on  the  main  trends  of  the  Ey  IXBz,  and  A//  traces.  It  is  worth 
mentioning  that  an  examination  of  the  high-resolution  raw 
data  (not  shown)  does  not  change  our  results.  We  consider 
first  dc  then  ac  aspects  of  these  observations. 
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In  the  interval  between  1 802  and  1 803  UT  A5z  increased  by 
about  50  nT,  indicating  that  was  about  0.1  %n  A/ni2. 
Simultaneous  variations  in  ABz  and  ABy  (Figure  2) 
support  our  adoption  of  the  infinite  current  sheet  approxi¬ 
mation  [Fung  and  Hoffman,  1992].  Combining  equation  (1) 
and  J^hm’s  Smiddy  et  al.  [1980]  obtained 
dy  yABz  —  j  =  0,  which  reduces  to 

ASz  —  1.256  •  E/>  •  =  const  (3) 

Equation  (3)  indicates  that  variations  in  ^z  and  ^  should 
track  each  otherjn  regiojig^of  constant  Up.  From  1802:10 
to  1802:35  UT  Ey  and  ABz  increased  by  ~35  mV/m  and 
~35  nT,  respectively.  During  thisjnterval  the  cross  corre¬ 
lation  coefficient  between  Ey  md  ABz  was  ~0.98.  Over  this 
interval,  one  can  estimate  Ep  as  ~0.8_rnhos.  Between 
1802:40  and  1803:10  the  slopes  of  the  ASz  and  Ey  traces 
have  opposite  signs  due  to  E/>  increase  to  1 .4  mhos.  At  the 
center  of  the  high-flow  channel,  where  «  55  mV/m,  the 

rate  of  column  Joule  heating  is  I^pE^^  «  2  mW/m^.  An 
integration  of  equation  (1)  along  the  spacecraft  trajectory 
shows  that  the  total  perturbation  is  proportional  to  J«  =  f  juds 
w  0.6  A/m.  "  " 

[17]  Figure  4  shows  perturbation  quantities  designated 
hEy,  bBz,  and  8A^  obtained  through  the  application  of 
bandpass  elliptic  filter  (0.05  <  /  <  0.45  Hz)  to  the  Ey 
ABz,  and  Nj  measurements,  respectively.  Spectral  analyses 
of  the  bEy  bBz,  and  bNj  measurements  are  shown  in  the 


Figure  2.  Top  and  second  panels  show  electron  and  ion 
energy-time  spectrograms  (logarithmic  scale);  1,  2,...  on 
the  color  bar  stand  for  10*,  10^. . .  l/(cm^  •  s  •  ster  •  eV). 
Third  panel  shows  horizontal  component  of  the  convection 
velocity  across  the  satellite  trajectory  in  meters  per  second. 
Bottom  panel  shows  deviation  of  measured  7(light  line)  and 
Z  (heavy  line)  components  of  magnetic  field  from  IGRF 
values.  The  dashed  line  indicates  the  SAPS’  poleward 
boundary. 


[is]  A  simple  application  of  Ampere’s  law  and  the 
current  continuity  equation,  using  the  infinite  current  sheet 
approximation  gives 


-  1  dABz 

/II  = 


Pp  dY  cos  ot 


=  V-U. 


(1) 


Here  po  is  the  pemieability  of  fi-ee  space,  Ij.  ~  Ep  AE  is  the 
height-integrated  ionospheric  current  perpendicular  to  Bp, 
and  Ep  is  the  height-integrated  Pederson  conductivity  (in 
the  iiightside  subaurqral  ionosphere  Hall  conductance  E^^  is 
negligible).  Herein,  4'  designates  spatially-averaged  (“dc”) 
perturbations  obtained  through  the  application  of  lowpass 
(At  >^20  s)  filter  to  ’^';yj|,  AS,  E,  and  Ep  are  expressed  in 
pA/m  ,  nT,  mV/m,  and  siemens,  respectively. 

[16]  For  a  satellite  traveling  poleward  at  a  speed  of  7.5 
km/s,  equation  (1)  is  approximated  by 


7|I«0.1-ASz. 


(2) 


UT 

Figure  3.  Variations  of  the  poleward  convection  electric 
field  (top),  deviation  of  the  eastward  component  of  the 
geomagnetic  field  from  IGRF  (middle),  and  ion  density 
(bottom). 
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Figure  4,  (top)  Oscillations  of  the  electric  (heavy  line)  and 
magnetic  (light  line)  fields,  (middle)  Ion  density  oscilla¬ 
tions.  (bottom)  Power  spectral  density  of  the  electric  field 
(heavy  line),  magnetic  field  (light  line),  and  ion  density 
(dotted  line)  oscillations. 


bottom  panel.  The  largest  spectral  component  for  the 
electric  and  magnetic  field  variations  is  at  <^0.17  Hz.  The 
bNi  spectral  peak  is  at  ~0.165  Hz. 

[18]  Using  measured  values  of  Ni  and  Bq  we  have 

calculated  the  local  Alfven  speed  Va  =  where 

M  is  the  ion  mass.  Current  voltage  sweeps  of  the  retarding 
potential  analyzer  on  FI 5  indicate  that  was  the  domi¬ 
nant  ion  species-  Within  the  low  density  region  of  the 
SAPS  structure  Va  «  1300  km/s.  In  the  top  panel  of 
Figure  5  we  compare  calculated  values  of  Va  with  those  of 
‘‘phase  velocities”  Vph  =  |6E'y/55z|.  In  regions  where  Vp^l 
Va  ^  associated  waves  are  almost  purely  electromag¬ 
netic  (likely  Alfven  waves).  Significant  deviation  from 
unity  indicates  the  presence  of  a  strong  electrostatic 
component  (note  that  points  with  Vp^lVA  >  10  were 
excluded), 

[19]  Finally,  the  bottom  panel  of  Figure  5  gives  calculat¬ 
ed  values  of  the  field-aligned  component  of  Poynting  flux 
iS^  =  hE'^Bzl\XQ  (in  [iW/m^)  carried  by  electromagnetic 
waves.  Positive  (negative)  values  indicate  that  wave-energy 
flows  into  (out  of)  the  ionosphere.  Apparently,  in  the  region 
of  highest  plasma  flow  near  tibe  center  of  the  structure,  wave- 
energy  flows  away  from  the  ionosphere  at  a  maximum  rate 
of  >^15  p,W/m^.  It  is  noteworthy  that  was  less  than 
0.15  •  hBz  during  this  period. 


[20]  An  examination  of  the  figures  yields  the  following 
four  observations: 

[21]  1.  During  a  ^50-s  period  after  ^1802:10  UT  the 
three  traces  show  very  similar  fluctuations, 

[22]  2.  The  Ey;  and  Nt  variations  are  anticorrelated  over 
the  interval. 

[23]  3.  Prior  to  1 802:20  the  phase  difference  between  hEy 
and  hBz  vary  from  ^>^0°  to  ^^^90°,  while  after  ^1802:42  UT 
hEy  and  hBz  vary  in  phase  with  each  other.  During  the 
intervening  period  they  are  !=5il80°  out  of  phase. 

[24]  4,  The  field-aligned  component  of  Poynting  flux  Sx 
is  directed  into  (out  of)  the  ionosphere  where  the  variations 
are  in  (out  of)  phase. 

[25]  We  have  surveyed  measurements  of  Vh  and  Ni  on 
DMSP  F12,  F13,  and  F14,  taken  in  the  evening  sector  within 
±30  min  of  the  SAPS  crossing  by  F15  shown  in  Figure  2.  In 
each  instance  we  found  signatures  of  the  SAPS  structure 
with  the  maximum  westward  flow  speed  near  52°  MLat, 
extending  for  at  least  3  hours  in  local  time.  Anticorrelated 
irregularity  structures  appeared  in  the  V^  and  Ni  traces. 
Spectral  analysis  of  these  data  showed  narrow  peaks  near 
^0.17  Hz.  Electromagnetic  interference  from  spacecraft 
contaminates  the  body-mounted  magnetometers  on  F12- 
F14  and  render  estimates  of  Poynting  flizxes  imreliable. 

[26]  Table  1  lists  the  universal  times  (UT),  magnetic 
latitude  (MLat),  and  magnetic  local  times  (MLT)  at  which 
maximum  westward  plasma  velocities  (Fmax)  were  observed 
by  FI 5  within  16  SAPS  stmctures.  Also  listed  in  Table  1  are 
the  corresponding  maximum  electric  fields  (E'max)j  the  fre¬ 
quency  of  irregularities  (/),  and  the  field-aligned  component 
of  Poynting  flux  (iS^)  deduced  from  electric  and  magnetic 
field  irregularities  within  SAPS. 

[27]  One  can  see  that  electromagnetic  substructures 
observed  by  FI 5  within  SAPS  events  were  detected  in  two 
parts  of  the  storm’s  recovery  phase  (Figure  1).  The  first 
group  (0400-0800  UT)  occurred  in  the  early  recovery  phase 
near  the  deepest  excursion  of  SYMH  and  the  second  rise  of 
AE  after  the  first  northward  turning  of  the  IMF  Bz.  They 
were  located  ^3°-5°  equatorward  in  magnetic  latitude  of 
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Table  1.  SAPS  and  Wave  Structures  From  DMSP  F15,  6 
November  2001 


0356  -46°  22.8  1.7  35  0.2  +1/-2.S 

0423  46°  20.7  1.5  '45  0.15&0.35  +1200/- 500 

0539  -47°  22.7  0.55  10 

0752  50°  20.5  1.45  50  0.16  +10/-1S 

1412  -53°  22.3  0.8  30  0.17  +1.8/-0.4 

1439  54°  21.0  0.5  15 

1554  -50°  21.8  1.4  45  0.12  +100/-20 

1621  51°  21.0  1.9  60  0.19  +1/-2 

1736  -52°  21.3  0.7  20 

1803  52°  213  1.75  55  0.17  -75 

1915  -53°  21.4  1.0  25 


UT.  Prior  to  0423:30  UT  both  bBz,  and  6£y  had  strong 
spectral  peaks  near  0.15  Hz.  While  ftNj  oscillations  occurred 
near  0.15  Hz,  the  bulk  of  its  spectral  power  was  in  the  0.22 
to  0.26  Hz  band.  After  0423:30  UT  the  spectral  range  of  hBz 
and  bNi  variations  broadened.  fiEy  spectral  power  decreased 
near  0.15  Hz  but  increased  dramatically  at  0.35  Hz.  The 
overall  effect  of  data  presented  in  Figure  7  suggests  that 
plasma  and  field  fluctuations  observed  prior  to  0423:30  UT 
are  similar  to  those  observed  during  the  other  SAPs  events 
(cf  Figure  4).  Variability  observed  after  this  time  were  of  a 
considerably  more  complex  nature. 

[3i]  As  with  the  SAPS  event  near  1800  UT  (Figure  5)  we 
have  compared  “phase  velocities”  obtained  from  the  ratios 


+0.5/-0.75 

+0.75/-1 

+13/-5 


the  second  (1400—2400  UT)  group.  Five  cases  with  no 
entries  in  the  last  two  columns  had  irregularity  amplitudes 
too  small  to  be  measured  by  DMSP  sensors.  These  were 
associated  with  values  of  Vtnax  1  km/s.  However,  in  one 
case  with  V^^ax  -  0.8  km/s  irregularities  were  detected.  In 
cases  with  detectable  irregularities  the  dominant  frequencies 
were  confined  to  a  fairly  narrow  band,  0.18  ±  0.03. 

[28]  The  final  column  of  Table  1  lists  the  field-aligned 
Po5aiting  fluxes  in  p.W/m^  calculated  using  simultaneous 
measurements  of  bBz  and  6£y  We  adopt  a  convention  using 
positive  (negative)  values  of  Sx  to  indicate  net  flows  of 
electromagnetic  energy  into  (out  of)  the  ionosphere.  The 
range  of  Sx  covers  more  than  three  orders  of  magnitude. 
Almost  half  of  the  SAPS  events  with  detectable  irregulari¬ 
ties  had  energy  flows  >10  pW/m^,  associated  with  > 
45  mV/m.  Thus  the  example  presented  above  is  typical  of 
the  discussed  phenomena. 

[29]  The  obvious  exception  is  the  SAPS  event  encoun¬ 
tered  near  0423  UT  in  which  the  peak  to  peak  variability  of 
Vff  oscillations  exceeded  dc  plasma  velocities  and  the 
downward  Poynting  flux  reached  1.2  mW/m^.  Figure  6 
shows  particle  fluxes  as  well  as  the  measured  electric  and 
magnetic  fields  and  ion  density.  The  vertical  dashed  line 
marks  the  low-latitude  boundary  of  auroral  electron  precipi¬ 
tation.  Attention  is  directed  to  three  features  of  the  data. 
First,  significant  fluxes  of  precipitating  30  keV  ions  were 
detected  between  0423:20  and  0423:30  UT,  equatorward  of 
the  auroral  bormdary.  These  fluxes  probably  represent  the 
low-energy  portion  of  the  ring-current  distribution  function 
[Smith  and  Hoffman,  1974].  Second,  F)/ modulations  appear 
to  change  at  ~04:23:30  UT.  Finally,  the  background  plasma 
density  decreased  from  2  10^  cm“^  at  0422:00  to  a 
minimum  of  ~5  lO**  cm"^  near  0423:00  UT  then  recovered 
at  0424:30  UT.  Densities  sampled  between  0423:00  and 
0424:00  UT  were  marked  by  irregularities  similar  to  those 
found  in  Figure  3.  The  average  Alfv6n  speed  in  the  trough 
decreased  from  about  900  to  700  km/s  between  0423  00  and 
0424:00  UT. 

[3o]  Figure  7  displays  spectral  analyses  of  hNi,  hBz,  and 
hEy  measurements  taken  between  0423:00  and  0424:00  UT. 
For  simplicity  we  divide  the  interval  up  into  half-minute 
intervals  before  (Figure  7a)  and  after  (Figure  7b)  0423:30 
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Figure  6.  Top  panel  shows  electron  and  ion  differential 
number  fluxes  in  (cm^  •  s  •  ster  •  eV)“'.  (heavy  line) 
and  (triangles)  are  the  fluxes  of  1  keV  electrons  and  30 
keV  ions,  respectively.  '  and  stand  for  the  averaged 
electron  fluxes  in  the  channels  1 1  to  19  (>1  keV)  and  1  to  9 
(<1  keV),  respectively.  Second  panel  shows  horizontal 
component  of  the  convection  velocity  across  the  satellite 
trajectory  in  meters  per  second.  Third  panel  shows  the 
deviation  of  the  Z  component  of  the  geomagnetic  field  from 
IGRF.  Bottom  panel  shows  the  ion  density  variation.  The 
dashed  line  indicates  the  plasma  sheet  equatorward  boundary. 


MISHIN  ET  AL.:  SUBAURORAL  WAVE  STRUCTURES 


SIA  1-7 


Figure  7.  Power  spectral  densities  of  (dots),  hBz  (light 
lines),  and  8E'y  (heavy  lines)  (a)  before  and  ^)  after  0423:30 
UT. 


of  hEy  and  hBz  with  Alf/en  velocities  computed  from 
measured  values  of  Ni  and  Bq.  These  comparisons  show  that 
in  significant  portions  of  the  0423  to  0424  UT  interval  the 
perturbation  were  Alvenic.  DMSP  measurements  suggest 
that  intense  electric  field  events  associated  with  SAPS 
structures  are  highly  localized  and/or  short  lived.  They  were 
not  encountered  while  crossing  the  conjugate  SAPS  struc¬ 
ture  at  0356  UT.  No  large  amplitude  8E'y perturbations  were 
detected  by  F12-F14  in  SAPS  structures  crossed  within  a 
half  hour  of  0423  UT.  The  event  however  is  not  unique.  We 
have  observed  similar  8£'yOScillations  during  other  magnetic 
storms  on  the  poleward  side  of  SAPS  structures.  Like  the 
event  shown  in  Figure  6,  they  too  were  associated  with  the 
precipitation  energetic  ions  just  equatorward  of  the  auroral 
electron  boundary. 

4.  Discussion 

[32]  The  previous  section  provides  two  detailed  examples 
of  wave-like  variations  in  plasma  density  as  well  as  mag¬ 
netic  and  electric  fields,  observed  by  DMSP  FI 5  while 
crossing  midlatitude  troughs  during  the  magnetic  storm  of 
6  November  2001.  In  both  cases,  relatively  weak  field- 


aligned  currents  into  the  ionosphere  are  associated  with 
intense  poleward  electric  fields  that  drive  westward  plasma 
convection  within  the  troughs.  Data  summarized  in  Table  1 
indicate  that  these  wave-like  perturbations  are  stable  fea¬ 
tures  in  storm  time  midlatitude  troughs  and  are  found  in  the 
conjugate  ionospheres.  Comparisons  of  hEy/bBz  ratios  with 
local  Altven  speeds  indicate  that  both  electrostatic  and 
electromagnetic  variations  were  sampled.  Similar  FACs, 
electric  fields,  and  density  troughs  were  detected  by  F14 
and  FI 3  near  2000  and  1800  LT,  respectively.  fNf  and 
hEy  variations  had  spectral  characteristics  similar  to  those 
observed  at  the  local  time  of  FI 5. 

[33]  The  wave-like  perturbations  appear  consistently  in 
the  evening  MLT  sector  just  equatorward  of  the  auroral 
electron  precipitation  boundary.  Within  this  region  DMSP 
detected  weak  field-aligned  currents  into  the  ionosphere, 
strong  poleward  electric  fields,  and  deep  troughs  in  the  local 
ionospheric  plasma  density.  Low  values  of  Up  inferred  firom 
the  ratio  of  quasi  dc  magnetic  to  electric  field  variations 
suggest  that  the  troughs  extend  in  altitude  to  the  peak  of  the 
F  layer.  The  absence  of  sunlight  and  measurable  fluxes  of 
auroral  particles  indicate  that  the  E  layer  contains  little  to  no 
plasma. ' 

[34]  We  distinguish  between  two  manifestations  of  wave 
phenomena.  The  first  and  more  common  appears  in  spectra 
accumulated  during  the  intervals  1802:00  to  1803:00  UT 
(Figure  4)  and  0423:00  to  0423:30  UT  (Figure  7a).  They  are 
characterized  by  electromagnetic  and  plasma  density  varia¬ 
tions  that  have  the  same  fi*equency  range  of  0.1  to  0.2  Hz. 
The  second  and  more  dramatic  manifestation  of  wave-like 
perturbations  was  sampled  in  the  30  s  after  0423:30  UT  and 
was  characterized  by  large-amplitude  plasma  and  field 
oscillations  over  a  broader  range  of  fi-equencies  extending 
across  the  0.2  to  0.3  Hz  range.  In  the  second  case  the 
perturbation  densities  and  fields  appear  to  have  different 
firequency  responses.  In  this  and  other  magnetic  storms 
strong  waves  are  associated  with  the  precipitation  of  ener¬ 
getic  ions. 

[35]  The  remainder  of  this  discussion  considers  several 
possible  causal  mechanisms  for  inteipreting  the  first  type  of 
the  wave  phenomena.  It  is  necessary  to  examine  the  nature 
of  the  electromagnetic  fluctuations.  We  show  that  the  data 
are  consistent  with  a  complex  coupling  between  the  iono¬ 
sphere  and  magnetosphere  consequent  to  the  formation  of 
deep  troughs  associated  with  SAPS  structures.  We  also 
suggest  that  the  source  of  the  density  irregularities  was 
located  far  to  the  east  of  the  DMSP  trajectory. 

4.1.  Cavity-Mode  Interpretation 

[36]  It  is  tempting  to  suggest  that  DMSP  observations  of 

hBz,  hEy  and  hNi  variations  with  similar  fi^equencies  across 
high-speed  convection  channels  in  the  trough  indicate  the 
presence  of  resonant  shear  Alfven  waves  in  a  cavity.  This 
interpretation  is  consistent  with  observing  similar  variability 
in  the  conjugate  ionospheres.  By  comparing  the  ratio  hEyihB 
with  the  local  Alfven  speed  we  established  that  the  observed 
waves  have  both  electromagnetic  and  electrostatic  compo¬ 
nents.  The  electromagnetic  (EM)  part  propagates  along 
magnetic  field  lines.  The  field-aligned  component  of  the 
wave  vector  k\\  =  varies  with  the  magnetic  field  and 

plasma  density  along  individual  magnetic  field  lines.  We 
view  the  electrostatic  (ES)  variations  as  zero-frequency 
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spatial  structures  with  wave  lengths  \es=  VsJfES-  For  = 
7.5  km/s  andTi^s = 0. 1 5  Hz,  = 50  km  and  /fcy  «  0. 1 3  km“  ‘ . 

[37]  Data  presented  in  Figure  3  indicate  that  the  plasma 
and  field  oscillations  occurred  between  1802.00  and 
1803^:00  UT.  During  this  interval  DMSP  F15  moved  from 
48.4°  to  52°  in  magnetic  latitude.  In  a  magnetic  dipole 
approximation,  this  corresponds  to  X  =  2.27  and  2.56, 
respectively,  where  X  =  1/cos  is  the  Mcllwain  parameter 
and  A  is  the  invariant  latitude.  It  is  easily  shown  that  the 
length  of  dipole  field  lines  between  two  conjugate  iono¬ 
spheric  points  is  D(A)  =  2LRe  dC, cos  CVl  +  3  sin^  C, 

where  ^  is  the  magnetic  latitude  and  Re  6380  km. 
Performing  the  integration  gives  D(48.4°)  =  26730  km 
and  D(52  )  =  32350  km.  To  satisfy  the  resonance  condition, 
■^(A)  must  be  equal  to  an  integral  number  of  half  wave¬ 
lengths  between  the  conjugate  ionospheres.  If  these  distan¬ 
ces  represent  a  full  wavelength,  then  with /em=  0.16  Hz,  the 
average  Alfv6n  speed  between  the  conjugate  ionospheres  is 
i;i(48.4°)  =  4275  and  Ki(52°)  =  5180  km/s.  These  are 
significantly  higher  than  Alfv6n  speeds  calculated  from 
plasma  densities  and  magnetic  fields  measured  by  DMSP 
in  the  topside  ionosphere  or  inferred  for  conjugate  plasma- 
spheric  parameters  [e.g.,  Chappell,  1972].  Thus  we  con¬ 
clude  that  the  observed  wave  phenomena  cannot  result  from 
a  simple  cavity-mode  resonance  of  a  ~0.15  Hz  wave. 

4.2.  Ring-Current  Source  Interpretation 

[38]  Kennell  and  Petschek  [1966]  demonstrated  that 
whenever  the  fluxes  of  energetic  particles  trapped  in  the 
inner  magnetosphere  exceed  certain  threshold  limits,  the 
amplitudes  of  background  electromagnetic  waves  grow, 
causing  the  gyroresonant  particles  to  diffuse  toward  the 
atmospheric  loss  cone.  Whistlers  and  Alfv6n  waves  are 
responsible  for  the  pitch-angle  diffusion  of  energetic  elec¬ 
trons  and  protons,  respectively.  The  presence  of  field- 
aligned  currents  across  the  region  of  interest  indicates  the 
presence  of  energetic  ring  current  ions  at  conjugate  loca¬ 
tions  in  the  magnetosphere  [Vasyliunas,  1970].  The  reso- 
.  nance  condition  is 


Here  M  and  M// represent  the  masses  of  ring-current  ions  and 
plasmaspheric  hydrogen  ions,  respectively.  Substituting 
g3n'ofrequencies  into  equation  (5)  shows  that  the  energy 
of  resonant  ring-current  ions,  independent  of  species,  must 
be  >1  MeV.  This  greatly  exceeds  the  several  tens  of 
keV  measured  ion  energies  of  storm  time  ring  current  ions 
near  their  earthward  boundary  [Smith  and  Hoffman,  1974], 
Thus  we  conclude  that  ring-current  particles  cannot  be  the 
direct  sources  of  the  AlfV6n  waves  characteristic  of  the  first 
type  wave  phenomena  (Figure  2).  This  conclusion  is 
consistent  with  the  observation  that  DMSP  did  not  detect 
fluxes  of  precipitating  1-33  keV  ions  during  such  events.  At 
X  2.5  precipitating  particles  should  have  equatorial  pitch 
angles  <12  .  At  the  location  of  DMSP  the  local  dip  angle  is 
about  67°,  and  the  atmospheric  loss  cone  has  a  width  of  about 
58°.  The  vertically  looking  SSJ4  detector  intercepts  ions  with 
pitch  angles  of  23  ±  2.5°,  well  within  the  local  loss  cone. 

[39]  It  is  worth  noting  that  lyemori  and  Hayashi  [1989] 
and  Erlandson  et  al.  [1993]  observed  down-going  electro¬ 
magnetic  waves  with  ~0.1-1  Hz  frequencies  (Pci  range)  at 
X  >  4  in  the  dusk  and  morning  sectors.  Erlandson  et  al. 
[1993]  suggested  that  those  are  electromagnetic  ion  cyclo¬ 
tron  (EMIC)  waves  generated  by  RC  protons  in  the  RC- 
plasmasphere  overlap  region.  Could  this  mechanism  con¬ 
tribute  to  formation  of  the  down-going  Poynting  flux  and 
energetic  ion  precipitations  in  the  dusk-evening  sector 
(Figure  6)?  Denton  et  al.  [1992]  have  shown  that  the 
fr®ju®ncy  of  the  most  unstable  EMIC  waves  varies  as 
Eeqinc  .  Given  Baj{L)  ~  X“^,  one  can  expect  at  X  ~  2 
the  frequency  of  order  several  Hz  for  EP  and  0.3  Hz 
for  O  .  Therefore  O  -driven  EMIC  waves  within  the  RC- 
plasmasphere  interaction  region  cannot  be  ruled  out  as  a 
possible  contributor  to  the  second  type  wave  phenomena 
(Figure  6  after  0423:30  UT). 


2x/-V||  =2iv/+ 

where  V||  and  i||  are  field-aligned  components  of  velocity  and 
wave  vector,  respectively,  and  is  the  local  ion  gyrofre- 
quency.  The  most  prolonged  interactions  occur  near  the 
magnetic  equator  where  the  plasma  density  and  magnetic 
field  slrengft  change  most  slowly.  In  the  magnetic  dipole 
approximation,  at  X = 2  near  the  magnetic  equator  B^q  w  3 800 
nT.  Energetic  ring  current  ions,  however,  may  be  H*  or  O* . 
Their  cyclotron  frequencies  in  a  3800  nT  magnetic  field  are 
57  and  3.6  Hz,  respectively.  As  for  a  0.16  Hz  wave /<  f", 
equation  (4)  can  only  be  satisfied  if  the  wave  and  energetic 
ion  travel  in  opposite  directions.  The  dispersion  relation  for 
Alfvdn  waves  is  ^i,  =  The  local  Alfvdn  speed  is 

determined  by  the  local  magnetic  field  and  the  density  of 
cold  plasmaspheric  ions.  Twical  plasmaspheric  densities 
near  X  =  2  are  ~  5  •  10®  m~^  and  are  mostly  ions 
[Chappell,  1972].  Substitution  into  equation  (4)  shows  that 
the  energy  of  resonant  ions  is 
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4.3.  Magnetosphere-Ionosphere 
Coupling  Interpretation 

[40]  Striations  of  relatively  high  and  low  plasma  density 
in  the  nightside  trough  mark  the  regions  of  correspondingly 
(4)  higher  and  lower  height-integrated  Pedersen  conductivity. 
To  maintain  current  continuity,  polarization  electric  fields 
must  develop  to  enhance/diminish  X^in  regions  of  negative/ 
positive  plasma  density  variations.  We  have  already  noted 
that  as  expected  variations  in  6Xy  anticorrelated  with  those 
of  hNi.  Altered  potential  distributions  near  hN,  variations 
must  then  map  to  the  magnetosphere.  The  role  of  Alfv6n 
waves  in  transferring  this  information  to  the  magnetosphere 
has  been  discussed  by  a  number  of  investigators  who  model 
auroral  arcs  as  latitudinally  narrow  strips  of  enhanced  iono¬ 
spheric  conductivity  [e.g.,  Ogawa  and  Sato,  1971;  Maltsev 
et  al,  1977].  In  its  proper  frame  of  reference,  equatorward 
convecting  plasma  experiences  time  varying  electric  fields 
(Vy  •  V)Ey.  Consequent  field-aligned  currents  are  carried  by 
Alfv6n  waves  that  distribute  space  charge  (V  ■ 
along  magnetic  field  lines  to  turn  them  into  the  appropriate 
new  equipotentials  [e.g.,  Kan  and  Sun,  1985].  These  models 
predict  that  as  rapidly  developing  bN/  structures  alter  poten¬ 
tial  patterns,  electromagnetic  energy  in  the  form  of  upward 
moving  Poynting  flux  should  be  observed  above  the  iono¬ 
sphere.  Consistent  with  this  prediction,  DMSP  detected 
(5)  episodes  of  upward  wave  Poynting  flux  while  crossing  the 
SAPS  stracture. 
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[4i]  Miura  and  Sato  [1980]  pointed  out  that  the  trans¬ 
mission  of  electromagnetic  energy  from  the  ionosphere  to 
the  magnetosphere  cannot  be  the  only  interaction.  Rather,  a 
feedback  develops  whereby  the  magnetosphere  also  affects 
the  ionosphere.  Altered  potential  distributions  in  the  mag¬ 
netosphere  affect  the  drift  velocities  of  the  energetic  particle 
populations  that  in  their  proper  reference  frames  see  time- 
varying  electric  fields.  The  consequent  inertial  drifts  carry 
currents  j_L  =  (MNie/B^)d8E/dt.  In  regions  where  the  ener¬ 
getic  particles  experience  dbWdX  <  0,  the  conditions  for  an 
MHD  generator  jj_  •  6E  <  0  are  satisfied.  That  is  such 
regions  extract  mechanical  energy  from  the  particles  to 
launch  Alfven  waves  carrying  field-aligned  current  and 
Poynting  flux  toward  the  ionosphere.  Qualitatively  this 
adaptation  of  the  Miura  and  Sato  [1980]  scenario  coincides 
with  the  DMSP  observations. 

4.4.  Source  for  Observed  Plasma  Irregularities 

[42]  The  DMSP  measurements  should  be  considered  on 
two  spatial  scales  sizes.  The  larger  scale  size  is  the  '^4.5° 
latitudinal  width  of  the  SAPS  structure  and  the  smaller  one 
is  that  of  the  plasma  irregularity  structures.  The  coincidental 
formation  of  plasma  density  troughs  and  SAPS  structures 
was  predicted  based  on  theoretical  considerations  by  Schunk 
et  al  [1976]  and  Banks  and  Yasuhara  [1978].  They  argued 
that  intense  electric  fields  in  the  subauroral  ionosphere 
would  lead  to  the  formation  of  deep  troughs  through  a 
combination  of  plasma  outflow  and  enhanced  recombina¬ 
tion  stemming  from  the  charge  exchange  reaction  +  N2 
— >  NO^  +  N  whose  rate  kjs/^  increases  in  the  convecting/ 
heated  plasma.  Actually,  Schunk  et  al  [1976]  made  use  of 
the  relation  k^^  ^  ^  [E  x  Bq]"*.  Anderson  et  al  [1991] 

have  argued  that  the  plasma  outflow  related  to  frictional  ion 
heating  is  the  main  cause  of  the  formation  of  the  deep 
troughs.  Moffett  et  al  [1998]  reported  enhancements  of  the 
ion  and  electron  temperatures  ^4000  and  6000  K  within  the 
SAPS-related  troughs.  It  is  known  that  increases  with 
the  level  of  vibrational  excitation  of  N2  (see  the  latest 
assessment  of  the  laboratory  measurements  by  Viggiano 
and  Williams  [2001]).  As  the  vibrational  population  of  N2 
rises  with  electron  heating  may  significantly  contribute 
to  the  trough  formation  [e.g.,  Moffett  et  al,  1998;  Pavlov  et 
al,  2000]. 

[43]  Importantly,  no  significant  (>100  m/s)  plasma  out¬ 
flows  were  detected  within  the  density  troughs  under  study. 
On  the  other  hand,  significant  variations  of  the  ion  and 
electron  temperature  were  observed  (Figure  8).  Note  that  the 
time  resolution  of  Tg  was  4  (top  panel)  and  32  s  (bottom), 
while  Ti  was  sampled  at  a  rate  of  0.25  s”^  One  can  see  that 
the  electron  temperature  achieved  maximum  9200  and 
6500  K,  respectively,  at  the  density  minimum.  Just  the 
opposite,  Ti  decreased  to  '^lOOO  K  at  the  center.  The 
electron  energy  flux  of  Qt  ~  (1-3)10^^  eVcm“^s“^ 
from  the  ring  current-plasmasphere  interaction  region 
[e.g.,  Moffett  et  al,  1998]  can  explain  the  electron  temper¬ 
ature  variation  Tg  oc  Qj/Nt,  In  turn,  the  ion  temperature 
decrease  in  the  center  can  be  understood  in  terms  of 
frictional  heating,  provided  the  plasma  density  in  the  heating 
region  decreases  stronger  than  ET^, 

[44]  Why  nearly  monochromatic  density  irregularities 
with  wavelengths  of  several  tens  of  km  should  form  is 
beyond  the  scope  of  this  observational  report.  The  likely 
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Figure  8.  Variations  of  the  electron  (dots)  and  ion 
(triangles)  temperature  withiu  the  (top)  0423  UT  and 
(bottom)  1803  UT  SAPS;  solid  line  shows  the  (scaled- 
down)  ion  density. 


source,  the  drift  instability,  requires  the  E  x  B  drift  to  have 
a  component  in  the  same  direction  as  the  plasma  density 
gradient  [e.g.,  Ossakow  and  Chaturvedi,  1979].  It  implies  a 
weak  (tangential)  westward  (eastward)  component  of  the 
electric  field  at  subauroral  latitudes  that  renders  the  equa- 
torward  (poleward)  wall  of  the  trough  drift-unstable.  To  the 
contrary,  the  most  intense  irregularities  were  colocated  with 
the  fastest  average  sunward  flows  in  the  ituddle  of  the^ 
trough.  One  can  speculate  that  the  reduction  of  die  plasma 
density  due  to  electron  heating  and  vice  versa  may  serve  as 
a  generation  mechanism  of  the  density  irregularities.  Indeed, 
the  relations  hNi  -hTg  with  u  >  0  and  7;  oc  QjINi  provide 
positive  feedback  in  the  system.  However,  a  mechanism  that 
defines  the  scale  size  of  the  thermal  perturbations  has  yet  to 
be  found. 

[45]  We  suggest  that  the  irregularities  observed  by  DMSP 
between  1802:00  and  1802:50  UT  were  not  created  locally. 
With  these  observations  in  mind,  it  is  useful  to  ask  where 
cold  plasmas  sampled  by  DMSP  F15  between  1801:20  and 
1802:50  were  located  1  hour  prior  to  observation.  To  answer 
this  question,  it  is  necessary  to  recall,  that  the  corotational 
motion  of  the  plasma  was  subtracted  from  SSIES  measure¬ 
ments  shown  in  Figure  2.  At  the  a  latitude  of  50°  the 
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corotation  speed  is  about  0.38  km/s.  Cold  plasma  detected 
before  (after)  1801:50  UT  reached  DMSP  F15  from  earlier 
(later)  local  times  to  the  west  (east)  of  its  trajectory.  In  a 
corotating  frme  of  reference  plasma  associated  with  the 
largest  amplitude  irregularities  moved  westward  about 
1.12  km/s.  Assuming  that  the  plasma  velocity  was  nearly 
constant,  at  1702:00  UT  the  irregularities  were  about 
4000  km  to  the  east.  At  50°  latitude  one  hour  in  local  time 
corresponds  to  an  east-west  separation  of  r\j  1 300  km.  Thus 
1  hour  prior  to  detection  the  plasma  element  would  be 
located  ~3  hours  to  the  east,  that  is,  in  the  midnight  sector. 

[46]  Empirical  models  of  plasma  convection  show  that 
plasma  streamlines  (equipotentials)  exit  the  polar  cap  cross 
the  auroral  oval  with  substantial  alignments  in  the  north- 
south  direction.  In  the  auroral  oval,  streamlines  develop 
larger  east-west  alignments.  Some  streamlines  continue 
toward  the  dayside  through  the  auroral  oval  while  other 
pass  to  subauroral  latitudes.  Near  1800  UT,  DMSP  F13 
(near  1800  MLT)  and  F15  (near  2100  MLT),  respectively, 
measured  potential  drops  of  about  -80  and  -70  kV  across’ 
the  afternoon  convection  cell.  At  both  locations  the  potential 
drop  across  the  SAPS  was  ~-25  kV.  Burke  et  al.  [1998] 
showed  that  during  magnetic  storms  the  potential  drops  at 
subauroral  latitudes  are  significant  fractions  of  the  total 
potential  in  the  afternoon  cell.  The  equipotentials  that  reach 
subauroral  latitudes  cross  the  polar  cap  near  the  separatrix 
between  the  afternoon  and  morning  convection  cells.  These 
considerations  suggest  the  following  scenario: 

[47]  1.  Enhanced  magnetic  merging  at  the  magnetopause 
captures  finite  volumes  of  dayside  plasma  then  convects 
them  as  plasma  patches”  into  and  across  the  polar  cap 
[Kelley  and  Vickrey,  1984]. 

[48]  2.  Within  the  polar  cap  the  trailing  edge  of  the 
patches  are  drift  unstable  yielding  alternating  regions  of 
relatively  high  and  low  density  at  F  layer  altitudes 
[Tsunoda,  1988;  Crowley,  1996].  Significantly,  Kivang  and 
Heelis  [1997]  reported  DE  2  observations  of  fully  structured 
patches. 

[49]  3.  After  convecting  across  the  polar  cap  the  striated 
topside  plasma  enters  the  nightside  auroral  oval  where  the 
plasma  density  is  augmented  at  E  region  altitudes  due  to 
auroral  precipitation.  The  main  effect  of  the  E  region  is  to 
enhance  transport  processes  in  the  F  region,  shorting  out 
polari^tion  electric  fields  [e.g.,  Farley,  1959],  suppressing 
the  drift  instability,  thus  freezing  existing  large-scale  stria- 
tions  in  place  and  eliminating  small  scales. 

[50]  4.  Plasma  drifting  along  equipotentials  near  the 
separatrix  between  the  dusk  and  dawn  cells  convects  to 
subauroral  latitudes  where  it  becomes  entrained  within  a 
SAPS  structure  in  the  midnight  sector.  In  the  absence  of 
sustaining  auroral  electron  precipitation,  E  layer  plasma 
quickly  recombines  and  disappears.  Within  the  SAPS, 
deep  troughs  are  formed  at  F  layer  altitudes.  In  the  process 
the  striations  that  developed  in  the  polar  cap  are  main¬ 
tained  as  alternating  streams  of  high  and  low  plasma 
density.  However,  since  the  directions  of  convection  and 
density  gradients  have  become  nearly  perpendicular,  no 
further  drift-mode  growth  occurs. 

[51]  5.  As  indicated  in  the  subsection  on  M-I  coupling,  the 
regions  of  alternating  high  and  low  electric  field  are  the  sites 
of  generation  of  Alfv6n  waves  carrying  field-aligned  cur¬ 
rents  to  the  magnetosphere  and  to  the  conjugate  ionosphere. 


[52]  It  is  easy  to  test  the  hypothesis  that  the  plasma  and 
field  striations  observed  by  F15  within  the  SAPS  originated 
in  the  polar  cap.  We  have  examined  plasma  densities 
measured  by  DMSP  F13  and  FI 5  in  the  polar  cap  near 
the  times  listed  in  Table  1.  During  periods  when  striations 
were  detected  in  the  midlatitude  trough  we  also  found 
stmctured  plasma  in  the  polar  cap.  Densities  within  the 
structures  were  ~(2.5  ±  1)  •  10^  cm~^.  If  our  conjecture  is 
correct  the  plasma  density  at  840  km  must  have  decreased 
by  a  factor  of  10  in  convecting  from  ~85°  MLat  in  the 
midnight  sector  of  the  polar  cap  to  r\j  52°  MLat,  2100  MLT 
in  the  subauroral  ionosphere. 

5.  Conclusion 

[53]  Oscillations  in  electric  and  magnetic  fields  and  plas¬ 
ma  density  were  observed  by  DMSP  satellites  in  the  evening 
sector  at  subauroral  latitudes  during  the  6  November  2001 
magnetic  storm.  Those  were  detected  within  fast  convection 
streams  near  the  deepest  depletions  of  the  midlatitude  trough. 
The  latter  coincide  with  the  electron  temperature  enhance¬ 
ments  >6000  K.  We  distinguish  between  two  manifestations 
of  wave  phenomena.  The  first  and  more  common  is  charac¬ 
terized  by  electromagnetic  and  plasma  density  variations 
within  the  same  frequency  range  of  ~0.15  Hz  in  the 
spacecraft  frame  of  reference.  The  second  type  is  character¬ 
ized  by  large  amplitude  plasma  and  field  oscillations  over  a 
broader  r^ge  of  frequencies  ~0.2  to  0.3  Hz.  The  perturba¬ 
tion  densities  and  fields  appear  to  have  different  fi-equency 
responses.  In  this  and  other  magnetic  storms,  strong  waves 
are  associated  with  the  precipitation  of  30  keV  ions. 

[54]  The  observed  waves  eomprise  electromagnetic  (in 
part  Alfv6nic)  and  electrostatic  modes.  Electromagnetic 
energy  fluxes  can  be  directed  both  into  and  out  of  the 
ionosphere.  The  density  perturbations  appear  to  be  extended 
east-west  corrugations  in  the  plasma  flow  streams  with 
north-south  wavelengths  of  ~50  km.  The  hEy  and  6Ai 
variations  were  anticorrelated,  as  required  for  current  con¬ 
servation.  Alfv6nic  perturbations  are  consistent  with 
CTpected  effects  of  irregular  potential  distribution  around 
ionospheric  density  irregularities  mapped  to  the  magneto¬ 
sphere.  Inertial  currents  act  to  generate  mesoscale  field- 
aligned  currents  carried  by  Alfv6n  waves,  as  was  previously 
discussed  with  regards  to  auroral  arcs  formation.  On  the 
other  hand,  O*  -driven  EMIC  waves  within  the  RC-plasma- 
sphei^e  interaction  region  cannot  be  raled  out  as  a  possible 
contributor  to  the  second  type  wave  phenomena.  We  sug¬ 
gest  that  ENi  irregularities  observed  by  DMSP  satellites  in 
the  evening  sector  began  as  striated  plasma  patches  in  the 
polar  cap  that  convect  to  subauroral  latitudes. 
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